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Abstract. We model in detail a flare observed on Proxima Centauri with the EPIC-PN on board XMM-Newton 
at high statistics and high time resolution and coverage. Time-dependent hydrodynamic loop modeling is used 
to describe the rise and peak of the light curve, and a large fraction of the decay, including its change of slope 
and a secondary maximum, over a duration of more than 2 hours. The light curve, the emission measure and 
the temperature derived from the data allow us to constrain the loop morphology and the heating function and 
to show that this flare can be described with two components: a major one triggered by an intense heat pulse 
injected in a single flaring loop with half-length ~ l.Ox 10 10 cm, the other one by less intense heat pulses released 
after about 1/2 hour since the first one in related loop systems, probably arcades, with the same half-length. The 
heat functions of the two loop systems appear be very similar: an intense pulse located at the loop footpoints 
' followed by a low gradual decay distributed in corona. The latter result and the similarity to at least one solar 

event (the Bastille Day flare in 2000) indicate that this pattern may be common to solar and stellar flares, in wide 
Q ■ generality. 
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1. Introduction enduring events are generally more gradual and longer- 

lasting than compact flares, but exceptions exist. 
Coronal flares are known to be very complex phenomena, 

! and to involve multiple coronal structures, multiple spec- The soft X " ra y l[ & ht curve of flares consists, generally, 

tral bands and multiple physical mechanisms at a time. 01 a stee P risin S P hase ' a well-defined peak and a slower 

Furthermore, it is very difficult to define a typical coronal " generally exponential - decay. A gradual rise or a decay 

flare pattern (e.g. Golub & Pasachoff 1997). A "standard" composed by segments with different e-folding times (e.g. 

classification, of solar coronal flares divides them into 0stc n & Brown 1999) can also occur. 

two main categories, based fundamentally on the topol- Stellar flares are spatially unresolved and we have no 

ogy of the involved structures: compact flares and long- direct information on the morphology of the coronal struc- 

enduring events (Pallavicini et al. 1977). Compact flares tures involved, except in the presence of eclipses during 

occur mostly inside single loops whose shape and volume the flare (Schmitt & Favata 1999). The similarity of solar 

do not change significantly during the flare. Long-enduring and stellar X-ray flares, however, suggests that also stellar 

events, instead, occur in loop arcades, and higher and flares involve plasma confined in closed structures. 

higher loops are typically involved as the flare progresses. Empirical methods have been developed to infer the 

The arcade footpoints, best seen in the Ha line, appear gize of the flaring structures from the e -folding decay time 

as two ribbons getting more distant with time. Long- of Ught curves (Kopp & p oletto 1984; white et al 19g6} 

Send offprint requests to: F. Reale Poletto et al. 1988, van den Oord & Mewe 1989, Pallavicini 

* Based on observations obtained with XMM-Newton, an et al. 1990, Hawley et al. 1995, Reale et al. 1997, Reale 

ESA science mission with instruments and contributions di- & Micela 1998, see Reale 2002 for an extensive review 

rectly funded by ESA Member states and the USA (NASA) of these methods) . In the hypothesis of flares occurring 
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inside closed coronal structures, the decay time of the X- 
ray emission roughly scales as the plasma cooling time. In 
turn, the cooling time scales with the length of the struc- 
ture which confines the plasma: the longer the decay, the 
larger is the structure (e.g. Haisch 1983). A loop thermo- 
dynamic decay time has been derived (van den Oord & 
Mewe 1989, Serio et al. 1991) as: 



T~th 



120L q 



(1) 



where Lg and T7 are the loop half-length and the max- 
imum temperature of the flaring plasma, in units of 10 9 
cm and 10 7 K, respectively. The timescale above is de- 
rived under the hypothesis of impulsive heat released at 
the beginning of the flare. However, a significant heat re- 
leased during the decay may increase the decay time, and 
therefore lead to overestimate the loop length, if not cor- 
rectly diagnosed (Reale et al. 1997, Reale 2002). By means 
of extensive hydrodynamic simulations of decaying flaring 
loops, Reale et al. (1997) derived an empirical formula 
for the loop length, combining the information from the 
light curve and the trajectory of the flare in the density- 
temperature diagram : 



La = 



120/(0 



/(C) > 1 



(2) 



where tlc is the decay time derived from the light curve. 
This formula can be obtained from the expression of the 
loop thermodynamic cooling time (Eq^l, but includes a 
non-dimensional correction factor /(C), larger than one 
(i.e. a shorter loop length) if significant heating is present 
during the decay. The slope ( of the decay path in the 
density-temperature diagram (Sylwester et al. 1993) is 
maximum (~ 2) if the heating is negligible in the decay 
and minimum (~ 0.5) - the slope of the loci of the hydro- 
static loops with decreasing temperature - if the heating 
dominates the decay. This approach has been tested on a 
sample of solar flares observed with Yohkoh/SXT (Reale 
et al. 1997) and extensively applied to flares observed on 
stars of various spectral type (Reale & Micela 1998, Favata 
et al. 2000, 2001, Maggio et al. 2000, Gudel et al. 2001). 

The empirical methods are of easy application and ap- 
propriate to infer the size of the flaring loops and some 
information on the flare heating, provided that the light 
curve and the temperature and emission measure diagnos- 
tics are available with enough photon statistics and time 
resolution and coverage to derive a decay trend. 

An XMM-Newton observation of the nearest star 
Proxima Centauri, of spectral type dMe, includes a very 
well-observed flare, already presented in Gudel et al. 
(2002, hereafter Paper I) and further analyzed in Gudel et 
al. (2003, hereafter Paper II). Other flares have been ob- 
served on Proxima Centauri and studied in detail (Haisch 



1 The square root of the emission measure can been used 
as proxy of the density to construct the density-temperature 
diagram. 



et al. 1983, Reale et al. 1988, Poletto et al. 1988, Byrne 
& McKay 1989). However, the large effective area and the 
high time coverage of XMM-Newton has allowed to collect 
data at an unprecedented level of detail which motivate a 
deeper analysis aimed at a higher level of diagnostics. 

Time-dependent hydrodynamic loop models have been 
shown to provide a good description of the evolution of the 
flaring plasma (e.g. Peres et al. 1987, Reale & Peres 1995, 
Hori et al. 1997). In particular, in the hypothesis of com- 
pact flares, it is customary to assume that plasma moves 
and transports energy only along the magnetic field lines, 
and to consider one-dimensional models (e.g. Nagai 1980, 
Peres et al. 1982, Doschek et al. 1983, Nagai & Emslie 
1984, Fisher et al. 1985, MacNeice 1986, Gan et al. 1991). 

The light curve of the flare observed with XMM- 
Newton shows a very peaked maximum and a globally 
slow decay, with changing e-folding time and even a well- 
defined smoother secondary peak. In this work, the flare 
is modelled in detail throughout the late phases, well after 
the second peak. The high quality of the data and their 
detailed comparison to the model results allow us not only 
to constrain the length of the main flaring loop, but also to 
diagnose the involvement of other structures in the flare, 
and to constrain the heating functions (intensity, temporal 
and spatial distribution) of all the flaring structures. 

In Sect. |21 we describe the observation and the con- 
straints on the modeling, in Sect. [3] we describe our mod- 
eling approach in detail; in Sect. the simulations per- 
formed and the results obtained are presented, in Sect. [5] 
the results are discussed and in Sect. conclusions are 
drawn. 

2. The observation 

The flare has been detected during the observation of 
Proxima Centauri made by XMM-Newton (Jansen et al. 
2001) on 2001 August 12, with a total exposure time of 65 
ks. Fig. 2] shows the flare light curve in the 0.15 - 10 keV 
band, collected, in small window mode, with the PN detec- 
tor (Striider et al. 2001) of the European Photon Imaging 
Camera (EPIC, Turner et al. 2001); MOS detectors are af- 
fected by pileup problems. High resolution X-ray spectra 
between 0.35 and 2.5 keV, taken simultaneously with the 
Reflection Grating Spectrometers (den Herder et al. 2001), 
are also available, and in particular OVII 22 A and Ne IX 
13.5 A He-like lines have been detected during the flare 
and analyzed in detail, yielding density estimates (Paper 
I). All data were analyzed using the XMM-Newton Science 
Analysis System (version 5.3, for RGS data version 5.4.1). 
The flare covers most of the final 20 ks of the observa- 
tion with a maximum luminosity £x,o.i5-io ~ 3.9 x 10 28 
erg/s. A large optical burst was captured with the Optical 
Monitor (Mason et al. 2001) in the rise phase of the X-ray 
flare, and may be a tracer of the production of non-thermal 
particles in the corona (Paper II). 

Fig. ^ shows the light curve of the first 10 ks of the 
flare, and the hardness ratio (ratio of 1-4.5 keV to 0.4- 
1 keV count rates) in the same time interval (see also 
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Fig. 1. Upper panel: Light curve (10 ks) of the flare on 
Proxima Centauri on 12 August 2001 as detected with 
the XMM-Newton EPIC-PN detector in the 0.15-10 keV 
band. The flare can be segmented into six phases, two 
rising (Rl, R2) and four decay ones (D1-D4), bounded by 
the vertical dashed lines. The solid lines mark the decay 
trends. Time t=0 corresponds to 17:00 UT of 12 August 
2001. Lower panel: hardness ratio (ratio of 1-4.5 keV to 
0.4-1 keV count rates) in the same time interval as the 
light curve. Time resolution is 300 s. 



Paper II). The count rate in the figure is the raw one 
extracted from the selected image region which includes 
90% of the total counts and should be further multiplied 
by the dead-time correction factor of 1.41. The bin size is 
20 s. The corresponding PN spectra, collected in 16 time 
intervals as shown in Fig. have been fitted with 2-T 
MEKAL models (Mewe et al. 1995) in XSPEC (Arnaud 



1996). Fig. |21 shows the values obtained for the dominant 
hotter component throughout the flare. 

The count rate reaches values as high as ~ 122 cts/s 
(after corrections). We can identify different phases of 
the light curve, which will be relevant for the modeling. 
During the initial rising phase (hereafter Rl), the emis- 
sion increases steeply by about two orders of magnitudes 
reaching a maximum in ~ 1 ks. The following decay is 
initially relatively rapid (Dl, with a duration of w 0.5 
ks) and then becomes slower (D2, w 1 ks). The time of 
the switch from Dl to D2 coincides with the time when 
the hardness ratio stalls and changes to a constant level 
(Fig. After « 2.5 ks since the beginning of the flare, 
the light curve rises again (R2), slowly, for 1.7 ks, reach- 
ing a second smoother peak at ~ 40 cts/s. The follow- 
ing decay is similar to the one after the first peak: fast 
first (D3, « 1 ks) and then slower (D4, « 5 ks). The 
light curve in phases Dl, D2, D3 and D4 can be reason- 
ably approximated with exponentials, with e-folding times 
Tm = 1.45±0.08 ks, t D2 = 2.66±0.12 ks t D3 = 1.62±0.06 
ks and tda = 4.35 ±0.06 ks, respectively, shown in Fig.^ 

Fig.|21shows that the maximum fit temperature of the 
hot component, \og(T b s ) « 7.4, is reached in phase Rl, 
somewhat earlier than the maximum of its emission mea- 
sure \og{EM) « 51.3. 

In the hypothesis that the bulk of the flare, the first 
peak, occurs inside a single flaring loop, the e-folding 
time of the light curve and the slope of the n-T path in 
the initial decay can be used to estimate the loop half- 
length, according to Eq. J5J), calibrated for the XMM- 
Newton EPIC/PN spectral response, already applied to 
a few events (Giidel et al. 2001, Stelzer et al. 2002, Briggs 
& Pye 2003), with: 

F(C) = C a 6Xp(-C/Ca) + Qa 

where 



(3) 



1.2 ±0.1 



c a = 11.6 ±0.5 Ca = 0.56 ±0.06 
0.4 < C < 1.9 

An expression for the loop maximum temperature 
T7 = T max /(10 7 K) can be derived from fitting hydro- 
static model loops with isothermal models: 

T max = 0.184 T 1 ' 15 " 



(4) 

We obtain \og(T max ) ps 7.6. 

From Fig. it can be noted that the initial decay Dl 
is very short and only two points are defined in the n-T 
diagram, too few to obtain a well-defined decay trend. If 
one assumes the initial decay Dl is entirely due to plasma 
cooling, with negligible heating, the empirical expression 
of loop length (Eqs. J5J, © and Q) yields an upper limit 
to the loop half length L up sa 1.6 x 10 10 cm. 

The emission measure shows a second peak during 
phase R2, as the light curve does, while the temperature 
is practically flat. Phase D3 includes just two points with 
large error bars. Phase D4 is better defined: the tempera- 
ture and EM of the hot component both decay monoton- 
ically, the temperature more slowly. 



4 



F. Reale et al.: Modeling an XMM-Newton Flare on Prox Cen 



7.6 



7.4 



7.2 



7.0 



6.8 



D2 f--4 
i-~t 

D3 



fR1 



6.6L 
24.8 



25.0 25.2 25.4 25.6 
0.5 log(EM[crrr 3 ]) 



7.6 
7.4 

K 72 

j? 7.0 

6.8 

6.6 
51.5 

E 
o 

5 50.5 

yi 

-5 50.0 
49.5 

































i « 1 


1, 






D1 


D2 


R2 


D3 


D4 : 



2000 



4000 6000 
Time (s) 



8000 



10000 



Fig. 2. Temperature and emission measure diagrams 
during the flare: the upper panel shows the density- 
temperature (n-T) diagram of the dominating component 
of the 2-T fitting, where EM 1 / 2 has been used as proxy 
for the density. The dashed line marks the evolution of 
the values and the end point of each phase is labelled. 
The lower panels show the evolution of the temperature 
and emission measure separately. Dashed lines as in Fig.Q] 



3. The modeling 

3.1. The set up 

The general approach to model this flare is an evolved ver- 
sion of the modeling of another flare observed on Proxima 
Centauri in 1980 with the Einstein Imaging Proportional 
Counter (Reale et al. 1988). At variance with the previ- 
ous modeling effort, here the modeling will include later 
phases of the flare and more than one loop component, 
allowing us to diagnose contributions of flaring structures 
other than the main loop and the heating function at late 
times. 



The model assumptions are those typical of a solar 
coronal flare loop modeling (e.g. Reale 2002): the flare in 
each loop is triggered by a strong heat pulse; the loop 
is initially at equilibrium (Serio et al. 1981) at the tem- 
perature (~ 4 MK) of an active region loop, not far 
from the peak of the EM distribution in quiescent con- 
ditions (w 3 MK in Paper II). The flaring plasma is de- 
scribed as a fluid confined in a closed semicircular loop 
with fixed geometry and constant cross-section, perpen- 
dicular to the stellar surface and unchanged during the 
flare. The plasma moves and transports energy only along 
the magnetic field lines running parallel to the loop, and 
can therefore be described with a single curvilinear co- 
ordinate. The plasma evolution is then described by the 
time-dependent hydrodynamic equations of mass, momen- 
tum and energy conservation as done in many previous 
works (see references in Section^), including, as significant 
physical effects, the gravity, the compressional viscosity, 
the radiative losses from optically thin plasma, and the 
thermal conduction. The stellar gravity and radius have 
been assumed g* = 10c/q and i?» = 0.15-Rq, respectively 
(Pettersen 1980, Segrensan et al. 2003). There are two 
external energy inputs: a low, constant and uniform one, 
which keeps the loop initially at equilibrium; a high and 
highly transient one, Q(s, £), which triggers the flare, and 
is assumed to be a separable function of space g(s) and 
time f(t) (e.g. Peres et al. 1987): 



Q(s,t) = H f(t) g(s) 



(5) 



where Hq is the peak value of the heating rate, s is the 
coordinate along the loop, t is the time. We consider 
Gaussian spatial distributions, centered on so and with 
width er: 



g(s) = exp 



so J 



2a 2 



(6) 



There is no reliable way to determine a priori the inten- 
sity, the spatial distribution, the duration and the time 
dependence of the heating function. We therefore proceed 
by educated guesses and refine the choices with the feed- 
back coming from the comparison of the data to the model 
results. We will consider, in particular, three alternative 
distributions g(s): two thin Gaussians centered at the foot- 
points, a single wide Gaussian centered at the apex, and 
a uniform heating [a L). As for the time dependence, 
we will consider a heat pulse, described as fit) = 1 for 
< t < 6tjj, and f(t) = at any other time. The heating 
decay is assumed exponential: 



f(t)=exp[-(t-t d )/T H ] 



(7) 



The time-dependent hydrodynamic equations have 
been solved using the revised version of the Palermo- 
Harvard numerical code with adaptive regridding (Betta 
et al. 1997, Betta et al. 2001). Symmetry with respect to 
the apex has been assumed, and a half-loop modelled. 
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A flaring loop model is set up by selecting the loop 
length and the heating function. The details and condi- 
tions of the loop before heat ignition are not critical for 
the simulation results, provided that the pressure is high 
enough to have a significant amount of mass in the chro- 
mosphere for evaporation (see Sect.^J. 

3.2. The analysis of the results 

The numerical solutions of the 1-D hydrodynamic plasma 
equations are in the form of plasma density, temperature 
and velocity distributions along the loop at progressing 
times. For comparison with observational data, from each 
density and temperature distribution, the plasma X-ray 
spectrum at the focal plane of the EPIC-PN detector is 
synthesized as done in several previous works (e.g. Reale et 
al. 1988, Reale et al. 1997, Reale & Micela 1998). We con- 
sider the MEKAL spectral code (Mewe et al. 1995) with 
a metallicity Z=0.5, as on average found in the spectral 
fits (Paper II). The MEKAL spectra are folded with the 
EPIC-PN response function used in Giidel et al. (2001). 
The final results are weakly dependent on the details (or 
minor changes) of the response function, since they are 
mainly based on the analysis of global observables, such 
as the light curve in a broad spectral band (0.15 - 10 keV). 
The normalization of the light curve obtained from the 
loop model to the observed one provides the loop cross- 
section area, which is a free parameter in the model. 

The loop model spectra are fitted with single temper- 
ature (1-T) model spectra (the same as those used to 
synthesize them). The fitting provides a best-fit "aver- 
age" temperature Tfn, and an analytical normalization 
factor, which, multiplied by the loop cross-section area, 
yields the emission measure. The 1-T fitting is performed 
in the 0.8-10 keV sub-band. This allows us to compare 
Tfn to the temperature of the hot component of the multi- 
temperature fitting of the data, at least for Tfu ^ 10 MK. 

Whenever fitting the observation data requires the 
combination of two model loops, we synthesize the total 
emission at a given time by summing the two focal plane 
spectra - one for each loop, with the appropriate cross- 
section area - at that time. We then analyze the resulting 
sequence of spectra, one for each time, as we do for single 
loop spectra: we derive the light curve and fit the spectra 
with single temperature models. 

4. The results 

The modeling of this flare will be described following the 
flare evolution. It will first address the flare peak, i.e. 
phases Rl and Dl in Fig. then the first decay (D2), 
and finally the second peak and the late decay (R2, D3 
and D4). 

4.1. The flare peak 

We model the initial and most intense phase of the flare 
with a single flaring loop; we will call it loop A. 



4.1.1. The length of loop A 

The modeling requires, first of all, that we set the loop 
length. As mentioned in Sectional the empirical scaling 
laws applied to phase Dl provide an upper limit L max = 
1.6 x 10 10 cm. In the lack of a well-defined path in the 
n-T diagram, and therefore of reliable information about 
the heating decay, any length shorter than this may be 
appropriate. 

We will show here results for three loop half-lengths, 
namely the upper limit L = 1.6 x 10 10 cm, an intermediate 
value L = 1.0 x 10 10 cm and the half-length obtained for 
the Einstein flare (0.7 x 10 10 cm, Reale et al. 1988). The 
initial base pressures are po = 3 dyne cm -2 for the first 
two, and po — 4.3 dyne cm~ 2 for the last length value. 

4.1.2. The heat pulse 

The flare peak is driven by a strong heat pulse. The time 
dependence of the heat pulse is described in Section 13.11 
The data indicate a very rapid increase of the tempera- 
ture and therefore an impulsive heating. Typically in flares 
(and in their simulations as well) the emission measure 
still increases well after the heating has been turned off 
(e.g. Svestka 1976). The temperature is a better tracer of 
the heating duration, because the efficient thermal con- 
duction makes it promptly decrease as the heating de- 
creases. Fig. |2 and the time evolution of the hardness ra- 
tio (Fig. ^) suggest a duration of the order of 500 s; the 
choice of a pulse duration Stu = 600 s is good for all our 
simulations of this flare phase. 

A hint for the pulse intensity comes from the flare max- 
imum temperature. It is reached after a few seconds and 
then remains steady as long as the heating is constant, 
because thermal conduction rapidly balances the heating. 
By applying the loop scaling laws (Rosner et al. 1978) with 
logT = 7.6 (Section^, we obtain that, if the heating were 
distributed uniformly in a loop of half-length 10 10 cm, its 
intensity would be of the order of: 

E H fa 10~ 6 T 3 - 5 Zr 2 w 4 erg cm" 3 s" 1 (8) 

For this phase of the flare, we have considered two 
alternative spatial distributions of the heat pulse along 
the loop: i) just above the footpoints (sq = 0.1L and 
a = 0.03L); ii) centered at the apex (so = L and 
a = 0.3L). The width of the heating distribution influ- 
ences the simulation results little. 

4.1.3. Modeling phases Rl and Dl 

We will not report on the whole exploration of the 
space of the model parameters that we performed, but 
only on some cases providing representative results. We 
will discuss results for the three loop lengths listed in 
Section 14.1.11 For the intermediate loop length (L = 10 10 
cm), we show results for two cases, i.e. a heat pulse con- 
centrated at the loop footpoints with a maximum intensity 
of H = 60 erg cm~ 3 s _1 , and a heating deposited at the 
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Fig. 3. Left: Fitting the observed light curve (data points) of the phases Rl and Dl with hydrodynamic simulations 
of a single flaring loop. The figure shows the light curves obtained from a loop with half-length f 10 cm heated at 
the footpoints (solid line) and at the apex (dotted line), from a loop with half-length 0.7 x 10 10 cm heated at the 
footpoints (dashed line), and from a loop with half-length 1.6 x 10 10 cm heated at the apex (dashed- dotted line). 
Right: corresponding paths in the EM 1 / 2 -T diagram, obtained from fitting the model and the observed spectra with 
isothermal models (only the first five data points are shown). 



loop apex with a maximum intensity Hq = 12 erg cm~ 3 
s _1 . For the smallest length, we show results for a heat 
pulse concentrated at the loop footpoints with a maxi- 
mum intensity of Hq = 85 erg cm -3 s~ 1 . For the longest 
loop, we show results for a heat pulse at the loop apex 
with a maximum intensity of Hq — 10 erg cm -3 s _1 . 

The evolution of the flaring plasma confined in a loop 
is well-known from extensive previous modeling studies 
(e.g. Peres et al. 1982). The global characteristics of the 
evolution do not depend on the details of the heating (see 
also Section 15. 5|) : the heat pulse makes the temperature 
increase up to several tens MK along the whole loop in a 
few seconds, due to the high plasma thermal conduction; 
the dense chromosphere at the loop footpoints is heated 
violently, and expands upwards with a strong evaporation 
front. The upcoming plasma fills up the loop, very dy- 
namically first and then more gradually, approaching a 
new hydrostatic equilibrium at a much higher pressure. 
The loop X-ray emission increases mostly following the 
increase of emission measure. As the heating stops (or 
decreases), the temperature promptly begins to decrease 
everywhere in the loop. As mentioned in Sect. 14.1.21 the 
emission measure peaks later and then decreases too, with 
the timescale of the plasma cooling. The different model- 
ing choices lead to different time scales, values of density 
and temperature, and to different details of the evolution, 
which, in turn, determine differences in the X-ray emission 
and its evolution. 

Fig-Elshows the light curves obtained from computing 
3000 s of plasma evolution for the four representative mod- 



els described above (two for the intermediate loop, one for 
the short loop and one for the long loop). The model re- 
sults are sampled with a minimum sampling time of 20 s. 
The model light curves are matched to the data by syn- 
chronizing the light curve maxima. The loop cross-section 
areas obtained from the normalization of the light curves 
are 3.6, 4.2, 1.1 and 4.0 xlO 18 cm 2 , corresponding to as- 
pect ratios R/L = 0.11,0.12,0.08 and 0.07, respectively, 
where R is the radius of the loop cross-section, assumed 
circular. 

The light curves all rise steeply during chromospheric 
evaporation, and the steepness decreases when the evapo- 
ration becomes more gradual. The maximum occurs about 
400 s after the heat pulse has stopped, and the decay fol- 
lows the decrease of the emission measure due to plasma 
cooling. 

The rising phase obtained from modeling the shorter 
loop 2 is too slow to fit reasonably well the observed one. 
The emission evolution obtained with the long loop is too 
gradual around the flare maximum, due to the longer time 
scales implied, and is not able to describe the sharp flare 
peak 3 . 

The light curves obtained with the intermediate loop 
fit better both the rising and the peak phase. The heating 
at the footpoints fits the maximum better than the heat- 
ing at the apex. The rise is slower with the shorter loop 
because the initial evaporation front takes less time to fill 

2 A similar light curve is obtained with heating at the apex. 

3 A similar light curve is obtained with heating at the foot- 
points. 
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the loop, and, since then, the density - and the X-ray 
emission - increases more gradually. 

The T vs EM 1 / 2 plot of Fig.EHshows the paths obtained 
from fitting the spectra of the four flare models, at various 
times, to isothermal model spectra, in comparison with 
fittings of the data (first five points). All models match 
reasonably well the first four data points, all included in 
phase Rl and Dl. They depart from the fifth data point, 
which belongs to a later phase, as the light curves do after 
time t ~ 1500 s, and exactly where the hardness ratio 
stops decaying and gets constant (Fig.QJ. 

The results shown so far suggest us that the model 
which best fits phases Rl and Dl of the flare is the one 
with the loop of intermediate length (L = 1.0 x 10 10 cm) 
and the heating at the footpoints. We will consider this as 
the starting point for fitting the following phases. 
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4.2. The decay phase D2 

In phase D2, the decay of the light curve slows down signif- 
icantly, as shown in Fig.^ This trend cannot be explained 
with the cooling of a longer loop (Eq. QJ), because the de- 
cay is initially faster. Nor can it be explained with a heat- 
ing gradually decaying from the peak value: there would 
be a single slower decay trend. 

This change of slope may be explained in two alter- 
native ways: i) in loop A, a low residual heating, much 
lower than the initial impulsive heating and with a differ- 
ent evolution, remains active; ii) another loop is beginning 
to flare and its rising emission overlaps the emission of loop 
A, slowing down the overall decay, and determining also 
the second flare peak (R2+D2). 

The latter hypothesis will be discussed in the next 
paragraph, together with later flare phases. We now ex- 
plore the hypothesis of the residual heating: it must be 
significantly less intense than the strong initial pulse, so 
as to become important only later in the decay, and must 
decrease gradually so to drive the late decay. We neglect 
the effect of such residual heating as long as the initial 
pulse is on. 

Fig. 01 shows the fitting obtained with three different 
decaying heating functions, each switched on at the end of 
the heating pulse deposited at the footpoints, i.e. ta = 600 
s and th — 4500 s in Eq. (JJJ: a) uniform in corona, and 
with intensity Hq — 1.2 erg cm~ 3 s~ x ; b) the same, with 
lower intensity Ho = 0.75 erg cm -3 s _1 ; c) at the foot- 
points, with the same spatial parameters as the impul- 
sive heating and with intensity H$ = 15 erg cm -3 s -1 . 
Integrating over the loop length, the total rates of the 
heating a) and c) are « 1/4 of the total rate of the im- 
pulsive heating, heating b) w 1/6. The loop cross-section 
area obtained to best-fit the light curve down to phase D2 
slightly changes to 3.3, 3.5, and 3.2 x 10 18 cm 2 , respec- 
tively. 

It is immediately apparent that the residual heating 
deposited at the footpoints cannot describe the light curve 
in phase D2. At time t ~ 2500 s a thermal instability oc- 



Fig. 4. Fitting the observed light curve (data points) of 
the flare phase D2 with a decaying heating switched on at 
the end of the heating pulse at the footpoints: uniform in 
corona and initial intensity 1.2 erg cm -3 s _1 (solid line), 
the same with initial intensity 0.75 erg cm -3 s" 1 (dotted 
line), and with a decaying heating deposited at the foot- 
points (15 erg cm -3 s _1 , dashed line). Data points as in 
Fig.0 



curs, and the light curve first increases and then suddenly 
drops: indeed any decaying heating deposited at the foot- 
points has been found to be unable to describe this decay 
phase, because a thermal instability invariably occurs. 

Phase D2 appears instead to be described more ade- 
quately with the residual heating deposited uniformly in 
corona and as low as 1 /4 of the impulsive heating rate at 
t = 600 s. An even lower heating rate makes the emission 
decrease too fast. An equal amount of heating more lo- 
calized anywhere in the coronal part of the loop (e.g. at 
the apex) does not bring significantly different results, and 
we will henceforth refer to a residual heating generically 
deposited in the corona. 

4.3. The second peak and late decay: phase R2, D3 
and D4 

After time t 2500 s, the light curve rises again to 
form the second lower maximum. The question, of course, 
is how this maximum is produced. Occam's razor argu- 
ment would suggest us simply that a second heating pulse, 
weaker than the first one, occurs in loop A around that 
time. However, we find that a second heat pulse as intense 
as to produce the necessary emission increase in the same 
loop invariably leads to a temporary temperature increase; 
such temperature increase is not observed (Fig.|2J). On the 
contrary, the temperature stays constant in phase D2 and 
slightly decreases in phase R2. A similar trend is present 
also in the time evolution of the hardness ratio (Fig. HJ. 
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Fig. 5. Fitting the observed flare light curve [data points) 
from phase Rl to D3 with a model consisting of the sum 
{solid line of a footpoint-heated flaring loop with half- 
length 10 10 cm (with no decaying heating, dotted line)), 
and of a second top-heated flaring loop with half-length 
2.5 x 10 10 cm, ignited at the same time as loop A (dashed 
line). The two lower panels show the emission measure and 
temperature versus time obtained from fitting the model 
and the observed spectra with isothermal models. Data 
points as in Figs.nand|21 

The only way to produce the second emission peak 
with no significant temperature change is that a second 
loop system gets involved in the flare, and its emission 
adds to the one of loop A. We call this second loop, loop B, 
and model it separately from loop A. In addition to its 
length and its heating function, for modeling loop B we 
have to set a time shift of the heating switch on with re- 
spect to loop A. It is not trivial to constrain all these 
parameters because the evolution of loop B must be de- 
coupled from the decay tail of loop A. 



We start noticing that phases D3 and D4 of the light 
curve (Fig. ^) are similar to phases Dl and D2: the decay 
is initially faster and then slows down. The e-folding time 
in phase D3 is only slightly longer than the one in the 
corresponding Dl phase (see Sect. [3 and Fig. 0}. We take 
this as an indication that loop B may be similar to loop A, 
i.e. same length, and we choose to check this assumption 
against the possibility of a much longer loop B, namely 
L=2.5 xlO 10 cm. Indeed, loop B could be even shorter 
than loop A, if a residual heating were present during 
phase Dl. However, this would imply two different regimes 
of residual heating, one in phase Dl and another in phase 
D2, and introduce another set of free parameters in the 
modeling, which we prefer not to do, if unnecessary. 

Since the slower evolution expected from a longer loop 
may naturally lead to a delayed emission peak, we have 
explored the possibility that the flare in this long loop 
B is triggered at the same time as that in loop A, with a 
lower intensity and longer duration. In this hypothesis, the 
slower decay D2 may be explained with the superposition 
of the continuation of the fast decay Dl with the rise of 
the flare of loop B (hypothesis (ii) in Sec. 14. 2|) . In this 
specific scenario, we will drop the decaying heating in loop 
A. 

FigEl shows the results obtained with a heating pulse 
located at the top of loop B with Ho = 1 erg cm~ 3 s -1 , 
constant for 3200 s. The figure shows the light curve of 
the two flare components separately, and the light curve 
obtained by summing the spectra of loop A and loop B 
(with distinct cross-section areas) at corresponding times, 
compared to the observed light curve. The light curve of 
loop B rises very gradually and peaks at time t ~ 3500 
s, more than 2000 s later than the flare in loop A. The 
loop cross-section of this second loop that best fits the 
light curve is 6.2 x 10 18 cm 2 , which corresponds to quite 
a small aspect ratio R/L w 0.06. Fitting the total spec- 
tra of loop A+loop B with isothermal models, the global 
evolution of the emission measure is well described. The 
evolution of the temperature instead shows a deep mini- 
mum at time t w 2000 s. This is not present in the data, 
and suggests us to reject this model. A similar tempera- 
ture dip appears also if such long loop B is heated at the 
footpoints. 

Fig. HJ1 shows results obtained with a loop B twin of 
loop A, and with a heating duration and spatial distribu- 
tion of loop B identical to that of loop A, but triggered 
2600 s later, with a rate ten times lower, Hq = 6 erg cm~ 3 
s _1 . In this alternative scenario, the decaying heating of 
loop A is maintained. The best-fitting cross-section area 
of loop B is 16.7 x 10 18 cm 2 , five times larger than the area 
of the first flaring loop, or, equivalently, an arcade of five 
loops equal to loop A. This combination describes better 
the temperature trend, but the light curve shows a small 
dip at time t ~ 3000 s. As suggested by the residuals, the 
dip can be filled - and the fit further improved - simply 
by adding a third minor flaring component, adjusting the 
loop cross-section areas and the heating time shifts appro- 
priately (Fig.EJ). The best combination that we find is to 
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Fig. 6. Fitting the observed flare light curve (data points) 
from phase Rl to D3 with a model consisting of the sum 
(solid line) of two flaring loop systems with the same half- 
length and similar heating function: one (loop A, dotted 
line) is heated 2600 s earlier and 10 times more intensely 
than the other (loop B, dashed line). A residual heating 
sustains the decay of both loops. Data points and lower 
panels as in Fig. 



add two loops B with area 9.0 x 10 18 cm 2 and 10.8 x 10 18 
cm 2 (2.5 and 3 times the area of loop A, still compatible 
with an arcade of ~ 5 loops equal to loop A), and heated 
with time shifts of 2200 s and 2800 s, respectively, since 
the start of the heating of loop A. 

In alternative, we may think to fill the gap of the light 
curve by considering a longer-lasting heating pulse in loop 
B, but we checked that this choice fails to reproduce ade- 
quately the temperature evolution. 

The latest decay D4 can be reasonably well fitted by 
assuming a residual heating of loop B (or of two loops B) 
with the same characteristics and e-folding time as the one 
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Fig. 7. Fitting the observed light curve (data points) of 
the whole flare including two other flaring loop compo- 
nents. Residuals (data counts minus model counts) are 
also shown. Data points as in Fig. 

of loop A, and the superposition of the decays of loops A 
and B. 



5. Discussion 

5.1. The general scenario 

This work describes the hydrodynamic loop modeling of 
an X-ray flare observed on Proxima Centauri with XMM- 
Newton. The data are very detailed: we can distinguish six 
well-defined phases in the light curve and a well-defined 
path in the density-temperature diagram. Our approach 
has been to model each phase in detail taking the time- 
resolved density/temperature information into account. 

We find that this flare is best-described with the fol- 
lowing combination of components: 

— The initial phase of the flare, including the rise phase, 
the main peak and the initial decay, occurs in a single 
loop, loop A, with half- length 1.0 x 10 10 cm. 

— This phase is triggered by a heat pulse deposited at 
the loop footpoints over ~ 10 min. 

— A (four-fold) weaker residual heating is left in loop A. 
It is deposited in the coronal section of the loop and 
decays slowly with an e-folding time of ~ 1 hour. 

— About half an hour after the first heat pulse has 
stopped, other minor heating pulses ignite other loops, 
probably an adjacent arcade, and produce a second mi- 
nor peak in the light curve. 

— The loop arcade is made of ~ 5 loops with same length 
and cross-section area as loop A. 

— The heating function of the arcade is very similar to 
that of the main flare. 
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Table 1. Parameters of best-fit loop models 



Parameters 


Loop A 


Loop B 


Geometry 


Half-length (10 lu cm) 


1.0 


1.0 


Cross-section (10 18 cm 2 ) 


3.3 (3.6) Q 


17 (20) b 


Aspect (R/L) 


0.10 (0.11) a 


0.22 (0.25)'' 


Morphology 


Single loop 


Arcade (~ 5 loops) 


Heat pulse 


Location 


footpoints 


footpoints 


Rate c (10 28 ) erg/s 


27 


14 


Start time (s) 





2600 (2200 - 2800) 6 


Duration (s) 


600 


600 


Total energy c (10 32 erg) 


1.6 


0.8 


Heat decay 


Location 


corona 


corona 


Initial rate c (10 28 ) erg/s 


7.2 


1.8 


Start time (s) 


600 


3200 (2800 - 3400) 6 


e-folding time (s) 


4500 


4500 


Total energy c (10 32 erg) 


3.2 


0.8 


Plasma parameters 


Max. temperature (MK) 


46 


20 


Max. apex density (10 11 cm -3 ) 


1.1 


0.2 


Max. velocity (km/s) 


1400 


800 



a - If no residual heating is included in the decay phase. 

6 - assuming that two arcades of loops B are ignited in a sequence (600 s one from the other). 
c - assuming a loop aspect 0.10 and 0.22, respectively. 



Table ^ summarizes the parameters of the model flar- 
ing loops. 

Hydrodynamic modeling allowed us to derive a very 
detailed scenario with qualitative and quantitative con- 
straints on the loop morphology and on the heating func- 
tion. The modeling may not be unique: we cannot exclude 
that other combinations of the parameters may be found 
with a significantly higher modeling effort, and a few pa- 
rameters are not totally constrained. 

On the other hand, we notice that such a complex and 
detailed event is reasonably well-described in terms of only 
two dominant loop components and a well-defined heat- 
ing function, valid for both flaring loops. This result may 
provide a general pattern for the interpretation of stellar 
flares, even those which show light curves more complex 
than a simple rise+decay. 

The complex light curve of this flare may be in part 
explained by the larger collecting area of XMM together 
with its capability of long uninterrupted observations with 
respect to previous satellites, and we may have missed it 
in other stellar flares just because of insufficient S/N ratio 
and time coverage. On the other hand, there are stellar 
flares observed with enough time resolution, coverage and 
statistics, which are less complex (e.g. van den Oord & 
Mewe 1989, Pallavicini et al. 1990), and also many solar 
ones (e.g. Sato et al. 2003). 

If the data quality were not so high, we would not have 
been able to distinguish so many details of the light curve 
and to address them one by one, and we would have lim- 



ited our analysis, for instance, to an overall application of 
the empirical scaling law (Eq. @, ©, 0) approximat- 
ing the decay to a single decay. We would have obtained 
a decay time r S i„ ~ 4.3 ks, and a slope in the n-T dia- 
gram C, S i n ~ 0.5. With log(T o6s ) w 7.4 (Fig.EJ, we would 
have obtained Lg sa 13, i.e. 30% larger than the best value 
derived with detailed hydrodynamic modeling. The agree- 
ment between the two approaches is relatively good, also 
considering that the slope Csm is close to the lower limit of 
the applicability of formula J3J , and therefore it is better 
to take Lg as an upper limit. 

5.2. The loop morphology 

The hydrodynamic evolution of the plasma confined in- 
side a single loop of total length 2 x 10 10 cm (loop A) is 
able to explain the initial phases of the flare. The later 
phases, and in particular the second peak, instead require 
the ignition of a second loop system. The modeling tells 
us that a second longer loop triggered simultaneously to 
loop A can fit the second peak, but not the slow monotonic 
temperature decay after the flare maximum. To fit both, 
it is necessary to assume a residual decaying heating in 
the coronal section of loop A, and an arcade of ~ 5 loops 
identical to loop A triggered ~ 40 min later. We come up 
therefore with a flare involving a system of almost identi- 
cal loops, a single one first, and an arcade later. 

We have also realized that there is at least one solar 
event which presents several analogies with this scenario: 
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the so-called Bastille Day flare (14 July 2000). This is quite 
an intense solar flare (GOES class X5.7) whose light curve 
in the A1.12 filter passband of the Soft X-ray Telescope on- 
board the satellite Yohkoh shows a clear bump after the 
main maximum (see Fig. 3 in Aschwanden & Alexander 
2001). The same figure clearly shows also that the bump 
is associated with the spectacular ignition of a long arcade, 
also detected by the TRACE telescope. All this may sug- 
gest a certain similarity of the loop morphology of this 
event with that on Proxima Centauri. Also the timing of 
the light curve phases is not tremendously different: the 
bump of the solar flare occurs about 600 s after the peak, 
the second maximum of the Proxima Centauri flare oc- 
curs 3000 s after the first one. The different delay may be 
linked to the scale size of the loops: the solar arcade loops 
are a factor ~ 4 shorter than the predicted stellar loops. 
We sketch a possible scenario of the flaring loop system 
on Proxima Centauri scaled to the resolved scenario of the 
solar Bastille Day flare in Fig. |H1 

The half-length of both loop A and loop B is found to 
be of the same order as the estimated radius of Proxima 
Centauri (L/R* ss 1). This length is neither very far (1.4 
times) from the length estimated for the flare observed 
with Einstein (Reale et al. 1988), nor very large in absolute 
value (the radius of Proxima Centauri is indeed a small 
one), and it is still reasonable as for surface coverage (as 
shown in Fig. [BJ. 

5.3. The flare heating 

The modeling has provided us with detailed information 
on the heating deposition, both as for the spatial distri- 
bution and for the temporal evolution. A heat pulse de- 
posited in the coronal part of loop A seems unable to fit 
the sharp peak of the light curve. A heating deposited at 
the loop footpoints is instead more successful. 

On the other hand a heating at the footpoints is unable 
to drive the observed slow late decay, which instead seems 
to require a coronal location. 

The presence of both a steep rise phase and a slow 
late decay therefore suggests that both kinds of heating 
depositions, one at the footpoints and the other in corona, 
must be at work. 

We can also infer the relative weight of the two heating 
components. The footpoint heating is more impulsive, i.e. 
intense and short-lasting (a few min). The other heating 
component is less intense (~ 1/4) and releases its energy 
over a much longer time scale (one hour). Such features 
seem to be traced also by the optical light curve (see Paper 
II) , which shows a sharp peak and a slower decay starting 
at ~ 1/4 of the maximum optical count rate. 

It is interesting to note that: a) the two components 
contribute to the flare with comparable amounts of en- 
ergy (2 — 3 x 10 32 erg); b) a similar combination of these 
two heating components (« 10 32 erg each) in the loop ar- 
cade (loop B) is able to explain the second flare maximum. 
Although we do not exclude that refining the heating func- 
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Fig. 9. Heating function of best-fit multi-loop model 
(shown in Fig. of the Prox Cen flare according to our 
modeling. The main loop A (black) is heated first (thick 
solid line) with a pulse at the footpoints, followed by a 
lower and gradual decay (thick dashed line) deposited in 
the coronal segment of the loop. The arcade of loops B 
(grey) is heated later with a heating function similar to 
that of loop A (thin solid and dashed lines). 



tion of the arcade may further improve the fitting of the 
data, we have shown that using simply the same time pa- 
rameter values of the heating function yields a satisfactory 
description of the flare. 

The global heating function which drives the flare evo- 
lution of Fig. is shown in Fig. Given the close simi- 
larity of the heating function of the two flaring structures, 
independent although probably adjacent, we may advance 
the hypothesis that the two heating components may often 
be both present in many flares, and may represent a gen- 
eral characteristics of solar and stellar flares. The idea that 
one heating mechanism is probably insufficient to explain 
coronal flares is not new (e.g. Peres et al. 1987, Masuda 
et al. 1994), but our analysis may provide a detailed and 
quantitative pattern to be explored in other events. 

Note, in particular, that the heating at the footpoints 
may be driven by high-energy electron beams precipitat- 
ing along the loop from a reconnection site high in the 
loop, as often mentioned in the literature (e.g. Masuda 
et al. 1994), and as also traced by the optical light curve 
(Paper II). 

5.4. The energy budget 

On the basis of the model which best describes the X-ray 
data on this flare, we can make some considerations about 
its energy budget. Considering the loop aspect Ra/L « 
0.1 for loop A and Rb/L rj 0.22 for loop B, we obtain that 
the maximum heating rates injected in the main loop and 
in the arcade are w 27 x 10 28 erg/s and w 14 x 10 28 erg/s, 
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Fig. 8. Sketch of the possible scenario of the flaring loop system on Proxima Centauri scaled to the Bastille Day flare 
on the Sun. The size of Proxima Centauri and the flare loops are on scale. 



respectively. For comparison, the maximum luminosity of 
the flare has been estimated to be L X fi. 15-10 ~ 3.9 x 10 28 
erg/s (Paper II), i.e. ~ 15% of the maximum heating rate. 
The rest of the heating rate goes into enthalpy, conduction 
of thermal energy downwards to the chromosphere, and 
radiation at lower energies. 

The rate of the uniform heating at the beginning of 
the decays are ss 7.2 x 10 28 crg/s and w 1.8 x 10 28 erg/s, 
respectively. If we integrate in time, we obtain that the 
total energy released over the decay phase (3.2 x 10 32 erg 
for loop A and 0.8 x 10 32 erg for loop B) is for loop A 
twice, and for loop B equal to, the energy released in the 
heating pulse. 

Summing over the whole flare we obtain a total flare 
thermal energy input of 6.5 x 10 32 erg, to be compared to 
a total energy radiated in the 0.15-10 keV of 1.5 x 10 32 erg 
(Paper I), i.e. ~ 25% of the total injected energy. For com- 
parison, the total energy (of the analyzed part) and the 



peak X-ray luminosity of the flare observed on Proxima 
Centauri with the Einstein satellite were ~2x 10 31 erg 
and w 1.2 x 10 28 erg/s, respectively, i.e. about 1/3 and 1/4 
of the total energy and peak rate of the heating used in 
the relevant hydrodynamic modeling (Reale et al. 1988). 
That flare appeared to be less energetic and therefore, con- 
sistently, softer and with less efficient thermal conduction 
than the flare analyzed here. 

For comparison with a solar flare, the total thermal en- 
ergy involved in the Bastille day flare has been estimated 
to be ~ 0.5 x 10 32 erg (Aschwanden & Alexander 2001), 
about 1/10 of the energy input of this Prox Cen flare. 

5.5. The plasma evolution and parameters 

The best fitting model provides us with an insight on the 
evolution of the plasma involved in the flare and con- 
fined in the flaring structures. Some features are shown 
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in Fig. EI] ( see a l so Table 1). The most dynamic and in- 
tense evolution of the flare occurs in loop A: the plasma 
temperature rises in few seconds from about 3 MK to 30 
MK, and to an absolute maximum close to 50 MK after 
the first 100 s; it then settles at about 40 MK during the 
time the heat pulse is on. After the first 10 s plasma be- 
gins to evaporate significantly upwards from the chromo- 
sphere, with velocities above 1000 km/s and making the 
coronal density increase by almost two orders of magni- 
tude to about 10 10 cm -3 . In spite of the high velocity, the 
relevant Doppler (blue) shifts would be difficult to be de- 
tected even in the case of a favourable loop orientation to 
the line of sight, because restricted to very hot lines (> 10 
MK) undetected with the RGS, and a very short time in- 
terval (the initial 2-3 min) of the flare, as typically occurs 
in solar flares (e.g. Antonucci et al. 1987). The blue-shifted 
line component is therefore lost, because highly diluted in 
the typical time intervals of photon integration 10 min 
(Paper I). 

After ~ 1 min the strong evaporation front reaches the 
top of the loop. Since then, the plasma continues to fill up 
the loop much more slowly, reaching a density above 10 11 
cm~ 3 in the corona. 

The loop system B undergoes a similar evolution, but 
significantly less dynamic; the maximum temperature is 
about 20 MK, the coronal density at its maximum is one 
order of magnitude less than in loop A, and the maximum 
velocity one half that of loop A. 

The modeling highlights the presence of very hot 
plasma components, with twice the temperature value 
than the one obtained from simple data fitting. 

This is further apparent from the total distribution 
of the emission measure versus temperature, EM(T), ob- 
tained by summing the contributions of loop A and arcade 
B, shown in Fig. El The distributions are averaged over 
time intervals corresponding to the ones of the distribu- 
tions obtained from data analysis, shown in Figs. 8 and 9 
of Paper II (intervals A to D). The EM(T) distributions 
of Fig. llll share global similarities with those derived from 
the data, in particular to those of Fig. 8 in Paper II: a dom- 
inant hot component (<; 30 MK) in interval A, a broader 
and cooler distribution in interval B, an even cooler distri- 
bution with a long cool tail in interval C, the appearance 
of a significant cool (~ 10 7 K) component in interval D. 
The latter cool component relates to the ignition of the 
loop arcade B. 

The differences between the distributions derived from 
the hydrodynamic modeling and those derived from the 
data are not surprising because the integral inversion tech- 
niques used to derive the distributions in Paper II are ill- 
posed. In spite of this, the comparison with hydrodynamic 
modeling clearly provides a key for the interpretation of 
the main features of the distributions obtained from the 
data. 

The agreement of the hydrodynamic modeling results 
to the data is further confirmed by the focal-plane EPIC- 
PN spectra synthesized from the model of Fig. for in- 
tervals A to D, compared to the observed ones (Fig. 1X2^1 . 
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Fig. 10. Evolution of the temperature, density and veloc- 
ity along the main flaring loop. Distributions along half of 
the loop are shown at the labelled times (s). The profiles 
of the second loop system at their maxima are also shown 
(dotted line). 



The general trends are well reproduced by the model spec- 
tra; discrepancies mainly concern the intensity of some 
line groups, mostly related to differences of metal abun- 
dances, which we assume to be Z=0.5 for all elements in 
the model spectra. The good agreement in the hard sec- 
tion of the spectra is a further proof of the presence of 
significant hot plasma components. 

There are no density-sensitive He-like triplets in the 
RGS band for very high temperature plasma, and there- 
fore it is difficult to diagnose the predicted density values 
of the hottest flaring plasma. Available density diagnostics 
for this flare comes from line ratios of O VII and Ne IX 
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Fig. 11. Emission measure distributions versus temper- 
ature (EM(T)) obtained from hydrodynamic modeling of 
the flaring loop system, shown in Fig. El The distributions 
are averaged over time intervals corresponding approxi- 
mately to those of the distributions obtained from data 
analysis in Paper II (intervals A to D). 
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Fig. 12. Spectra (solid line) synthesized from results of 
the best-fit hydrodynamic loop model (Fig. 0) integrated 
over time intervals corresponding approximately to those 
of the spectra shown in Paper II (intervals A to D). Data 
are overplotted for comparison. 

groups obtained with the RGS (Paper I) . The line analysis 
provides density values «4x 10 11 cm~ 3 with a large error 
bar for O VII, and between 10 11 and 2 x 10 12 cm~ 3 , for 
Ne IX, in a time interval around the flare peak. In approx- 
imately the same time interval, the loop model yields an 



average density of i=s 2 x 10 12 cm~ 3 , and w 10 12 cm~ 3 at 
the temperatures of 2 and 4 MK of maximum formation of 
the respective ions. For the Ne IX, the average density ob- 
tained from the model is compatible to the large interval 
allowed by the data, although the Ne IX derived densities 
are highly uncertain due to severe line blending (Paper II) . 
The values obtained from the modeling at 2 MK are quite 
higher than those derived from the analysis of the O VII 
line. This may be motivated as follows: the O VII lines are 
intensely emitted both from the flaring plasma and from 
the remaining quiet corona. What we detect is therefore 
the sum of the two contributions, and the density an aver- 
age of the 2 MK plasma of the flare and of the whole Prox 
Cen corona. If we assume an average density value of the 
quiet corona at 2 MK (e.g. ~ 5 x 10 9 cm~ 3 , compatible 
with the pre- flare density value shown in Paper I), we can 
infer the relative weight of the components contributing 
to form the average density value and derive an estimate 
of the plasma volume involved by the quiet corona com- 
ponent. Just to give an idea of the volumes involved, this 
component could be contained in a shell surrounding the 
whole Prox Cen star of thickness ~ 10 10 cm, i.e. of the 
order of the stellar radius. 

From the properties of the confined plasma, we can 
also infer some properties of the magnetic field around 
the flaring structures. The pressure of the flaring plasma 
confined in loop A reaches values of the order of 1000 
dyne cm~ 2 . In analogy with the derivation in Maggio et al. 
(2000), we find that, to keep plasma confined, a minimum 
magnetic field of ~ 150 G is required; in order to extract a 
total energy of ~ 6.5 x 10 32 erg in a volume of the order of 
5 x 10 29 cm 3 and to maintain confinement, the initial non- 
potential magnetic field should have been at least ~ 230 
G. 

6. Final remarks 

In this work we model a flare observed at a high level of 
detail with XMM-Newton on Proxima Centauri. The good 
time coverage and resolution and the high count statistics 
of the data has allowed us to obtain very detailed diag- 
nostics by means of specific time-dependent hydrodynamic 
loop modeling. The modeling has allowed us, on the one 
hand to synthesize in detail a wealth of observables for 
comparison with data, and on the other hand to obtain a 
deep physical insight in the evolution and distribution of 
the confined plasma. 

On the other hand, the constraints provided by the 
data allowed us to discriminate among different model 
choices, such as different loop lengths, the presence of more 
than one flaring loop, the location of the heat pulse and 
of the residual heating, their intensity, relative timing and 
timescales. This work indicates that both heating compo- 
nents are necessary ingredients to explain this flare, and 
that a second loop system, probably an arcade, is required 
to explain the secondary maximum. 

In spite of the high degree of detail and of the many 
distinct trends present in the flare, relatively few model 
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components were sufficient to match the data reasonably: 
a loop and an arcade of loops, all with the same length, 
ignited with some delay by an intense heat pulse at the 
footpoints and a gradual residual coronal heating. Indeed 
a solar flare having a similar evolution, the Bastille Day 
flare, indicates that the scenario of involved loops that we 
find is realistic. It has been recently shown that sections 
of the Bastille Day flare are better described with a model 
of several concentric loops than with a single-loop model 
(Reeves and Warren 2002). Our results suggest the funda- 
mental ingredients which govern the X-ray flare evolution 
arc the plasma confinement, a few dominant loop systems 
with a fixed length, and a well-defined heating function. 
Changes of magnetic topology within each loop system 
seem to have a small influence on the X-ray evolution, 
probably because most of them are limited to a relatively 
small fraction of the life of the flaring arcade (see Fig. 5 
in Reeves and Warren 2002). 

The loop morphology and the heating function show 
a well-defined pattern which may be applied to interpret 
other stellar flares. The heating pattern may be applied, 
for instance, to stellar flares with multi-slope decay (Reale 
2002 and references therein) and with secondary maxima 
(e.g. Poletto et al. 1988, Pallavicini et al. 1990). 

It will be certainly interesting also to revisit solar flares 
showing similar features under such perspectives and to 
explore the theoretical implications concerning, in partic- 
ular, the flare heating mechanisms. 
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